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1 Introduction 



The one-dimensional (ID) Hubbard model has been one of the most fundamental and favorite integrable 
models in the non-perturbative quantum field theory. It exhibits on-site Coulomb interaction and 
correlated hopping which will possiblely reveal a promising role in understanding the mystery of the 
high-Tc superconductivity. Since Lieb and Wu 0, in 1968, solved the ID Hubbard model with periodic 
boundary conditions (BC) by coordinate Bethe ansatz, there has been a great deal of papers devoted to 
the study of the model. A remarkable step was done by Shastry who showed that the Hamiltonian 
of the ID Hubbard periodic chain commutes with a one-parameter family of transfer matrix of an 
equivalent coupled symmetric XY spin chain and who also gave a direct proof of the integrability of 
the model by presenting the quantum i?-matrix. Later on, Wadati and coworkers j^, further studied 
its integrability in terms of quantum inverse scattering method (QISM) Very recently, Martins 

and Ramos proposed a desirable way to solve the eigenvalue problem of the transfer matrix of the 
ID Hubbard model with periodic BC by means of algebraic Beth ansatz. Their approach provides us 
a unified way to solve a wide class of Hubbard-like models ^] by algebraic Bethe ansatz. 

On the other hand, in recent years, there has been much interest in the study of the quantum inte- 
grable systems with open BC, i.e. the systems on finite interval with independent boundary conditions 
on each end. Due to the presence of the boundary fields which lead to a pure back-scattering on each 
end of the quantum chain and the exhibition of the quantum group symmetry by special choice of the 
boundary parameters make the system possess rich physical properties [11, 24] in thermodynamical 
point of view. A systematic approach to handle the open BC for ID integrable quantum chains was 
proposed by Sklyanin ||l^. A further extension of Sklyanin's formalism to deal with more general 



class of models associated with Lie (super) algebras was proposed by Mezincescu and Nepomechie [13]. 
We also remark that the coordinate Bethe ansatz for ID Hubbard model with integrable boundary 
conditions was studied in ]17]. By now, though there are several authors ]|lj, |l^, 17] have stud- 
ied the open BC for the ID Hubbard model, the algebraic Bethe ansatz solution have not yet been 
achieved. Actually, the diagonalization of the transfer matrix which provide us with the spectrum of all 
conserved charges should be more essential in studying finite temperature properties of the integrable 
models ]18, ^ than diagonalization of the underlying Hamiltonian. But as we know the reflection 
equations for the ID Hubbard model are much more involved and the quantum i?-matrix does not 
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have the additive property that make it difficult to built up the necessary commutation rules among 
the diagonal fields and creation fields. In this paper, we intend to generalize Sklyanin's formalism to 
solve the ID Hubbard model with open BC. The eigenvalue of the transfer matrix and Bethe ansatz 
equations for the model will be given. It will be found that the model exhibits a hidden XXX spin 
open chain which play a crucial role to solve the model. 

This paper is organized as follows. In section 2, we shall recall the main results about open BC 
for the ID Hubbard model in order to introduce the notations which shall be used in this paper. In 
section 3, we perform the algebraic Bethe ansatz approach for the model. In section 4, we formulate 
the nested algebraic Bethe ansatz for the hidden XXX quantum spin open chain and present our main 
results. Section 5 is devoted to the conclusion. 



2 The ID Hubbard model with boundary fields 



Let us consider the ID Hubbard model with boundary fields determined by the Hamiltonian [|Tj, 15, 16| 

N-l N , 



1. 



j=l s j=l 

+p+(2nit - 1) + p_(2nii " 1) + 9+(2?^iVT - 1) + ^-(2njvi - !)• 



(1) 



Here p± and q± are the free boundary parameters characterizing the boundary fields. The coupling 
U describes the on-site Coulomb interaction and Oj^ and ajg are creation and annihilation operators 
with spins {s =| or |) at site j satisfying the anti-commutation relations 

{ajs, ay^/} = {a]^, a],^,} = 0, {ajs, a],^,} = 5^/5,,', (2) 

and Tijs = CL^js(^js is the density operator. The Lax operator is given 1^, 0] by 



Cj{u) 



(3) 



where 



fjs = sin u — (sin u — i cos u)njs, gjs = cosu — (cos u + i sin u)njs 
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With the grading P{1) = -P(4) = 0, ^(2) = P{3) = 1 and the constraint condition 

sinh2/i(n) _ U 
sin 2u 4 ' 

the Lax operator (^) satisfies the graded Yang-Baxter algebra 



(4) 



ni2{u, v) T (n) T {v) =T {v) T {u)ni2{u, v), (5) 
where the monodromy matrix T(u) is defined by 

Tiu) = CN{u)---Ci{u), (6) 

and 

Tiu)=T{u)0sI;Tiu) = l0sT{u). (7) 
here (8>s is the super direct product: 

For our convenience in practical calculation, we display the associated quantum 7^i2(ti, f )-matrix in 



Appendix. One may show that TZi2{u,v) enjoys the following graded reflection equations |16| : 



ni2iu,v) K_ (u)7^2l(^;,-n) K_ (v) =K_ {v)ni2iu,-v) K_ (u)7^2l(-^' , -n), (8) 



7^; 



StlSt2 / 



21 (v, u) Kf^ (n)7^^^l^*^ (-n, v) Kf' (v) 
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>StiSt2 



1 



Sti 



>StiSt2 



(9) 



which ensure the integrability of the model (|l|) provided that 

fKl±{u) \ 

K2±{u) 

K3±{u) 

V KA±{u) J 



K±{u) 



(10) 



where p+ = P- = '?-/2 
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= A„ 






sinn 


K2^{u) 


= A_ 




sinu)(e-''(") cosn-e'^(")C- 


sinn 


K3^{u) 


= A_ 




sinu)(e~^(") cosn-e'^(")^_ 


sinn 


A"4_(u) 


= 




sinu)(e-^(") cosn + e'^(")^_ 


sinn 


= W2 











Kl+{u) = A+(e-'^(")^+cosn + e''(")sinn)(e'^(")^+cosn + e-''(")sinn), 

K2+{u) = A+ (e'*^") e+ cos n + e-'*^") sin n)(e'^(")e+ cos n-e-'*^") sinn), 

i^3+(n) = A+ (e''^^) e+ cos n + e-'^^") sin n)(e'^(")e+ cos n-e-''^^) sinn), (12) 

K4+{u) = A+(e''(")e+cosn-e-'*(")sinn)(e-^(")C+cosn-e'*(") 



sm n . 



Here X± and are arbitary constants describing boundary effects. St^ stands for the inverse of the 
supertransposition in the space a. The supertransposition is defined by 



We would hke to remark that Zhou first time gave a class of boundary i^±-matrices equivalent to 
(PH) and ( p^ ) in terms of QISM. Consequently, using Lax pair formulation, the author |jl^ presented 
two class of boundary i^±-matrices leading to four possible boundary terms in the ID Hubbard open 
chain Hamiltonian. While, Shiroishi and Wadati studied the open BC for the model in terms of 
the graded version of QISM and also presented two class of the solutions to the graded RE. The second 
solution to the graded RE and (^) permits the boundary fields with p+ = — p_ and g+ = —q- 
corresponding to magnetic boundary fields (see ref. |jl5|, |l^] ). In this paper, we restrict to study 



the chemical boundary fields (11) and (12) basing on the consideration that this kind of boundary 
conditions will bring us a simple boundary i^-matrix for the hidden XXX open chain. To other kinds 
of boundary conditions, of course, we may treat them in a similar way. It is found that the Hamiltonian 
(^ is related to the double-row monodromy matrix 

T(n) = StvoK+{u)T{u)K^{u)T-^{-u) (13) 



in the following way 

r(n) = ciu + C2U^ + cz{H + const.) -\- ■ ■ ■ (14) 

where q, i = 1, • ■ ■ 4, are some scalar functions of boundary parameters. Strg denotes the supertrace 
carried out in auxiliary space vq. 



|0). 



(15) 



3 Algebraic Bethe ansatz approach 

According to the algebraic Bethe ansatz, let us first choose the standard ferromagnetic pseudovacuum 
state |0)i 

/I \ / I 

as a highest vector, which corresponds to the doubly occupied state. Following the notation introduced 

in Q, we define the monodromy matix Tin) as 

/ Biu) Bi{u) B2{u) F{u) \ 

Ci{u) Aii{u) Ai2{u) Ei{u) 

C2{u) A2i{u) ^22 (n) E2{u) 

\C^{u) C4(n) Csjn) D{u) j 
( B (n) Bi {u) B2 iu) F (n) \ 

Ci (n) An (n) A12 {u) Ei {u) 

C2 {u) A21 (n) A22 {u) E2 [u) 

\C3iu) Ciiu) C5iu) b{u)) 



Tiu) 



(16) 



r-\-u) 



(17) 
(18) 



and 



(19) 



T{u)K-{u)T~^{-u) 

( B{u) Bi{u) B2{u) F{u) \ 

Ci{u) Aii{u) Ai2{u) Ei{u) 

(72 (n) A2i{u) A22{u) E2{u) 

\C;{u) C^{u) C^{u) b{u) J 

It is not difficult to show that T-.{u) also satisfies the RE (^). Acting T{u) and T~^{—u) on the 

pseudovacuum state 

|0) = ®ili|0)i, (20) 
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we have the fohowing properties (upon a common factor): 



B{u)\0) = B (n)|0) = ^ ^e^M") |o), 



smu 



D{u)\0) = ^(n)l0) = |^e2'^wj''|0), 

L cos u J 

Aaa{u)\0) = Aaa (n)|0) = |0), a = 1, 2, 

Aab{u)\0) = Aab {u)\0) = 0, a ^ b = 1,2, 

Baiu)\0) / 0,Baiu)\0)^0,a = l,2, (21) 

Ea{u)\0) / 0, Ea {u)\0) ^ 0, a = 1,2, 

F{u)\0) / 0,F{u)\0)^0, 

Ci{u)\0) = Ci{u)\0) =0,i = l,---,5. 

Using the properties (|2l|), and the Yang-Baxter algebra 

T (-n)7^l2(^i, -u) T {u) =T (n)7^l2(^^, -u) T {-u), (22) 
and after some algebra, one can obtain 

B{u)\Q) = W^{u)B{u) B {u)\0), (23) 

Aaa{u)\0) = {!^^^^^l^Biu)Biu)+W^{u)Aaa{u)Aaa{u)}\0), (24) 

pi(u, -u) 

D{um = l^—(K2^{u)-^^p^)j2 AM Aaa in) 
[Pa{u,-u) Pi{u,-u) 

_^ P3i'^,-u) ~ ^ Ty-(n)Z)(n) D (u) 1 |0), (25) 
pi{u,-u) J 

Baiu)\0) / 0, Eaiu)\0) ^0, a = 1,2, (26) 

Aabiu)\0) = 0, a/6 = 1,2, F(u) /O, (27) 

(5,(1^)10) = 0, i = l,---,5, (28) 

where 



W,-{u) = 1, (29) 
_ (e-2M«) +e2^(")) Sinn cos n(^_e''(") cos -u-e-^^") Sinn) ^^^^ 
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^ cos 2ti(e~^'^(") cos^ u — e^'*^") sin^ u) 

^g-/i(u)^_ cQg^ _ g/iH sinu)(e'^*^")^_ cosn + e~'^^'^^ sinn) ^^^^ 

We also notice that the operators Ba{u), Ea{u) and F{u) are stiU creation fields, otherwise, C'j(n) are 
the annihilation fields. Via the transformations 



i:,(n) = Aaa{u)-^^^^^B{u), (32) 
pi(n, -n) 

D'iu) = Diu)-^^^^B{u)-—^j2AUu), (33) 
pi(n, -n) P4{u,-u)^^ 



we may express the transfer matrix (|T3|) in the following way 
r(n) = StroK+ (n)T_(n) 

2 

= H^+(n)^(n) + H^2+WE^-W + ^/W^'(^), (34) 



a=l 



where 



I , , (e + e ( )) smncosnsm2n ^, , 

PFr(n) = , , ^ rr-^^; 5— /(n) X 

cos 2n(e^"'^"-' cos^ n — e sin n) 

sinn - e'^^") cosn)(e'^(")C+ sinn - e'^'^") 



cosn 



(e '^(")^_|_ cosn + e'^^") sinn) (e''(")i^+ cosn + e '*(")sinn) 
+ _ (e-2feW +e^fcW) sinn cosn 
^2 N - (e-2h(n) cos2 n - e2M«) sin^ n) ^ ^""^ 

(e'i-H^^ cosn — e~^^^^ sinn)(e'^*^")^+ sinn — e~'*(") cosn) 
(■g-fe(u)^^ cosn + e'^^") sinn)(e'^(")^+ cosn + e~'*(") sinn) 



with 



(35) 



(36) 



^+ _ (e-feWg^ cosn - e^^") sinn)(e^(")g+ cosn - e'^^") sinn) ^^^^ ^^^^ 
4 (e~'^(")i^+ cosn + e'^(") sinn)(e'^(")i^+ cosn + e~''(") sinn) 



/(n) = e-2^'^(") cos^^ n sin^^ n i^l_(n)E:i+(n). (38) 



Now we proceed the key step to built up the necessary commutation relations between the diagonal 
fields and the creation fields respectively. From the RE (|8|) and definition (|Tg|), after many steps of 
substitution, we can get the following important commutation relations: 

\^ f \ pi(n,n)pio(n, -n) - / , . /gn\ 

B{u)Ba{v) = — -Ba{v)B(u) +u.t., (39) 

Pi{v,-u)pio{-u,-v) 
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D {u)Ba{v) 

Ba{u) Bk{v) 



where 



with the weights 

a{u, —v) 
a (— —u) 



P7{u,-v)pg{-v,-u) ~ . ^ I + fAn\ 

r :—Ba[v)D (u) + U.t., (40) 

pg{u,-V)ps(U,v) 

"^^^'°/^r^^(n, -v) {-V, -u)B,{v)X^,{u) + u.t., (41) 
(u,v) " ^ 



Pl{u, -v)pg 

pioju, -v)pi{-v, -u) 
pi{u,v)pio{v,-u) 
Pq{u,-v 



{b,{v) Bd{u) 



pio{u,-v] 

P6{V,-U) 

pio{v,-u) 

P8iv,-u)pQ{-V,-u) 
Pl0{v,-u)p8{-V,-u) 



■F{v) ( {I(^A{u))j . r {-v,-u) 
F{u) e (/0i(f)) 

[F{v)Biu)-F{u)B{v)].t. 



r{u, —v) 



r {—V, —u) 



C= (0,1,-1,0), 
/I 0\ 

a(n, —v) b{u, —v) 

6(n, —v) a{u, —v) 

Vo 1/ 

/I 0\ 

a{—v,—u) b{—v,—u) 

b {—v,—u) a{—v,—u) 

VO ly 



P3{u,-v)p4,iu,-v) - 1 



pg{u,-v)pio{u,-v) 

Pbi-v, -u)p8{-v, -u) + Pei-v, -u) 



b{u, —v) = 1 — a{u, —v), 



(42) 



(43) 



(44) 



b {-V, -u) = 1- a {-V, -u). 



(45) 



(46) 



Pa{-v, -u)ps{-v, -u) 

In the commutation relations (|39|)-(p|), we had to omit all unwanted terms because they take a big 
space to display. It turns out that the auxiliary matrices r{u, —v) and r {—v, —u) are nothing but 
the rational i?-matrices of isotropic six-vertex model. The structure of the auxiliary matrics is very 
important to solve the Hubbard-like |l^, |2l| models with open BC that exhibit a similar structure of 
the auxiliary matrices like Eq. (O) and (^3). If performing the parameterization introduced in f^]. 



one may also find 



U u + V 

a{u,-v) = ^ , b{u,-v) = ^ (48) 

u + v + u u + V + U 

a{-v,-u) = bi-v,-u) = - —- (49) 

u — V + u u — V + U 



In view of the commutation relation (^), the creation operators Ba,Ea do not interwine. So it is 
reasonable that the eigenvectors of the transfer matrices are generated only by the creation operators 
Ba{u) and F{u) or Ea{u) and F[u). Unfortmiately, it seems to be very difficult to construct the explicit 
form of the multi-particle vector even in the case of the Hubbard periodic chain |^. But it does have a 
similar recursive relation as that for the Hubbard periodic chain. Here we prefer the n-particle vector 
in a formal form, namely 

I ^n{vi,---,Vn)) = ^n{vu---,Vn)F''^'- -'''-\^). (50) 

Where the n-particle vector $„(7;i, • • • , f„) may be given by a recursive relation 

^n{vi,- ■ ■ ,Vn) = Ba{vi) ®^n~l{v2.- ■ ■ ,Vn) 

n _^ 

)B{vj)gf\{vi,---,Vn) (51) 

i=2 

- ®F{vi)]^n-2{v2, ■ ■ ■ ,Vj-i,Vj+i, ■ • •,«„)(/ (g) A{u))h'f\{vi, ■ ■ ■ ,Vn). 

i=2 

From the commutation relation (^2|), we can conclude that $„(ui, • • • ,f„) also satisfies the symmetry 
relation 



^n\Vi,- ■ ■ ,Vj,Vj+i,- ■ ■ ,Vn) — 7 r z 

Pl{Vj,Vj+i)pio{Vj+i, -Vj) 

^n{vi, ■ ■ ■ ,Vj+i,Vj, ■ ■ ■ ,Vn). r {-Vj+i, -Vj) (52) 

basing on the following relation: 



Pi{-Vj+1:-Vj) J _ - ^ Ps{Vj+l, -Vj)p6{-Vj+1, -Vj)p8{-Vj, -Vj+l) 

pi{vj,Vj+i) ■ J+ ' J p8{-Vj+i,Vj)ps{vj,-Vj+i)pe{-Vj,-Vj+i) 



This symmetry giving a restriction to the functions h^"^i{vi, ■ ■ ■ ,Vn) and gj^i{vi, ■ ■ ■ ,Vn) is very useful 
to deduce the coefficients and simplify the unwanted terms in the eigenvalue of the transfer matrix. In 
fact, after performing three-particle scattering, the explicit form of these coefficients can be fixed. But 



checking three-particle scattering is indeed a extremely tough work, we had to leave the coefficients to 
be determined. Explicitly, we display the two-particle state 

^2{VI,V2) = BaM)®Ba,{v2)- ~''^\ f {v^) ^ {I ® A{v2)) 

Pl0(f2, -Vl) 

Pl0{V2, -Vl)Ps[-V2, -Vl) 

In above expressions, F""^'' ■ ' > are the coefficients of arbitrary linear combination of the vectors reflect 
the " spin" degrees of freedom with = 1,2 and ^ plays the role of forbidding two spin up or two spin 
down at same site. F(u) creats a local hole pair with opposite spins. Acting the diagonal fields on the 
vector dl^), we way phenomenologically get 

B{u)\^,,{vi,---,vn)) = B{u)f[ P^i^^^^)P^o{u,-v.) |^^(^^^...^^^))^^.^.^ (55) 
D{n)\^r.ivi,---,vn)) = D'{u) n - ^'^""l -^^)Pf-J^i^ I ,^^(^^, . . . ^^^)) + u.t., (56) 

fj{ P9[U,-Vi)p8[U,Vi) 
a' f \ \ ( W a' t \T\ P4:{-Vi,-u)pio{u,-Vi) ^ 

ruiu - viY^i'^'ruiu + V2)lfg\ri2iu - V2)T,'l^ ' ' ' 

ruiu + Vn)tZl-'n2{u - Vn)Zl:i„ I ^n(t^l, " " " , Vn)) + U.t. (57) 

It follows that 

f \ \ ff. r w \Ti/ \TJ Pi(.^i^'^)Pio{u, -Vi) 

t{u) \ ^n{vi,- ■ ■ ,Vn)) = <W^{u)B{u)[[ 



pi{vi, -u)pio{-u, -Vi) 



provided that 



_lw+/ \n'r \TJ p7{u,-Vi)pg{-Vi,-u) 
+WZ{u)D [u] [[ — — 

P9{u,-Vi)ps{u,Vi) 

f=l Pl{u,-Vi)pg{u,Vi) J 

\^n{vi,--- ,Vn)) (58) 



'^"^"^#^U^.= A«(n = ^.{i).}). (59) 



i = 1, ■ ■ ■ n. 

Here ri2{u) = P.r{u) and A*^^^(n, {vi}) is the eigenvalue of the nested transfer matrix (|6l]), i.e. 

T'^^\u,{vi})F^^"'^" = AW(u,{i}i})F^i---^", (60) 
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where 

r(^)(n, {v^}) = TtoT^^\u)T^^^~\-u). (61) 
The nested monodromy matrices T^^\u) and r*^-*^^ ^(""w) read 

rW(^) = r^,{u + n)l\l-..ru{n + Vr.)l:';-\ (62) 
rW(-S) = ri2(n-^i)^:;::^,„...n2(n-tii)^^^ (63) 

So far, the eigenvalue problem of the ID Hubbard model with boundaries reduces to solve the nested 
auxiliary transfer matrix (6C) which corresponds to an isotropic six- vertex model with open boundary 



conditions. 

4 The nested Bethe ansatz 

In this section, we proceed the diagonalization of the auxiliary transfer matrix (|6l|). Following 



Sklyanin's formalism [12|, performing the nested Bethe ansatz have not been a difficult problem yet. 
It is easy to check that the ri2(ti)-matrix satisfies the Yang-Baxter algebra 

12 2 1 

ruiui - U2) T(^) (ni, {i,}) T(^) (us, =T^'^ (^2, {^i}) T^'^ {ui,m)ri2{ui - U2), (64) 

and the reflection equations 

1 22 1 

'^i2('Wi - U2) K^^ (ni)ri2(ni -|- U2) K^-^ (""2) =K^y {u2)ri2iui + U2) (ni)ri2(ni - U2), (65) 



ri2iu2 - ui) K^^^ {ui)ri2{-ui - U2 - 2U) {U2) 

2 1 
--K^^^ {U2)ri2{-ui -U2- 2U) {ui)n2{u2 - ui). (66) 



For our case, the {u) = I. Let us define the nested monodromy matrix 

r(^)(n)=r«(n)T«-'(-^)= (67) 



which also satisfies the RE ([65|). Using the main ingredients (p^)-([67|) describing the open BC com- 
patible with the integrability of the model and following all steps solving XX Z open chain in [l^, one 
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can present the following results: 

M 



provided that 

TT (ttj + U){uj -Vj + U) _ yj (uj + UI + 2U){uj -ui + U) 

~J[ {Uj + Vi){Uj -Vi) {Uj + Ul){Uj -ui-U) 

1 = 1, 

J = 1, • • • M, 

which indeed ensure the cancellation of all unwanted terms in (|6^). Here the "spin" part of the 
multi-particle states is given by 

I {v,})) = • • • S«(nAf)|0)(i), (70) 

where M is the number of holes with spin down, n is the total number of the holes. 
Finally, if we adopt the variables z±{vi) used in 0,i.e. 

^ COST^^2hK)^ ^^(^^) = !^e2M..), (71) 

sm Vi COS Vi 

and make a shift uj = Xj — U/2, the eigenvalue of the transfer matrix (|13|) is given as 

r(n) I ^n{vi, ■ ■ ■ ,Vn)) = 

+ W, [u)W, [u)[z^{u)\ 11 ^^^^^^^ _ _ ,^^u)/z^{vi)) 

. , 2(n + ^) sin^ n(l + + 1/ z^{vi)z+{u)) 

+ [u)W^ [u) 2-^^ l\cos^u{l-z^{v,)/z^{n)){l-l/z^{v,)z^{u))'^ 

^ {u + \i-U/2){u-\i-U/2) 
}X{u-\i + U/2){u + \i + U/2) 

_LW+r ^ 2-u 1^ sin^^(l + + 

+ H/2 i^)2n + C/lJcos2n(l-z_K)z+(n))(l-z+(tx)/z-K)) 

,y (n-A, + C//2)(n + A, + [//2) j' ^ ^ 
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provided that 



CK,C+)Cw,C-)[2-K)] =11 



.1 {vi-Xi + U/2){di + Xi + U/2) 



n 



CXj + Vi + U/2)CXj-Vi + U/2) 



\{Xj+v^-U/2){Xj-Vi-U/2) 



M 

n 



1=1, 

1+3 

j = !,■■■ M, i=l,--- ,n, 



CXj + Xi + U)CXj-Xi + U) 
CXj + ~Xi-U)CXj -Xi-U)' 



(73) 
(74) 



where 



e ^(")^_|_ sin u — e^^"^^ cos u 



± COS u — e' 



-h{u) 



smti 



(75) 



If we express the variable Z-{ui) in terms of the momenta ki (hole) by z^{ui) = e^^\ from the relation 
(p^), the energy is given by 



E^ = ^_+^^- {N/2 - n)U - ^ 4cos h. 



(76) 



i=l 



Equations (72)-(76) constitute our main results of this paper. It is found that the boundary fields are 
indeed nontrivial to the ground state properties and the boundary energy of the model. 



5 Conclusion 

We have formulated the algebraic Bethe ansatz solution for the ID Hubbard model with open bound- 
aries. Bethe ansatz equations, the eigenvalue of the transfer matrix and energy spectrum have also 



been given. Comparing our results with coordinate Bethe ansatz solution |17], the Bethe ansatz equa- 



tions (|73|) and (74) coincide with ones obtained in [0], apart from this, we presented explicitly the 
eigenvalue of the transfer matrix and conjectured the main structure of the n-particle eigenvectors. The 
results obtained provides us with a start point to study the thermodynamical properties of the model 
||6|, Especially, the proposed way is available to formulate the algebraic Bethe ansatz for other ex- 
tented Hubbard models with open BC, such as ID Bariev open chain p6| , 10|, Uq[Osp{2\2)] electronic 
system [20, etc. We also notice that if we add the chemical potential term i^X^j^i X]s('^js ~ 1/2) 
to the Hamiltonian (|l|), the integrability of the model require the associated quantum i?-matrix |22], 



which does not have crossing unitarity, should satisfy new RE. But the new class of the boundary 
ir-t-matrices |23] shall not change the Bethe ansatz equations (^) and (^). Nevertheless, if we add 
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the Kondo impurities |24, 25| JJ2ss' '^I'^ss'^s' ^° each boundaries, the model is also integrable with 
a certain boundary ii'-t -matrices which lead to new Bethe ansatz equations. We hope following this 
paper we shall present a class of integrable Kondo impurities for the ID Hubbard model in near future. 
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Appendix 



We display tiie quantum lZ{u, u)-matrix of the ID Hubbard model below ^ 
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with the Boltzmann weights 



pi = (cosn cosv e +sinu sinn e~ )p2, 
P4 = (cosn cos t> e~' + sin t; sinu e'')p2, 
Pq = (sinn cos ne~' — sin u cosn e')/92, 

pio = (sinn cosn — sinn cosn e~^)p2, 

(cosn cosn — sinn sinn e~') 

P3 = 5 ~ P2, 

cos^ u — sm n 

(cosn cosn e~' — sinn sinn e') 

P5 = 5 ~ P2, 

cos^ n — sm n 

e~''(cosn sinn — sinn cosn e~') 

P6 = 5 — P2, 

cos^ n — sm v 

and 

P8 = Pi - Ps; P7 = P4 - P5, 

I = h{u) - h{v), h = h{u) + h{v) 



15 



which enjoy the following identities: 

PAPl + P9PW = 1, PlP5 + P3P4 = 2, 
pI = P3P5 -'i-, pI = P9P10 + P7P8- 
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